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ABSTRACT 

We have made an asteroseismic analysis of the variable blue stragglers in the open cluster 
M67. The data set consists of photometric time series from eight sites using nine 0.6-2.1 me- 
ter telescopes with a time baseline of 43 days. In two stars, EW Cnc and EX Cnc, we detect 
the highest number of frequencies (41 and 26) detected in 5 Scuti stars belonging to a stellar 
cluster, and EW Cnc has the second highest number of frequencies detected in any S Scuti star. 
We have computed a grid of pulsation models that take the effects of rotation into account. 
The distribution of observed and theoretical frequencies show that in a wide frequency range 
a significant fraction of the radial and non-radial low-degree modes are excited to detectable 
amplitudes. Despite the large number of observed frequencies we cannot constrain the funda- 
mental parameters of the stars. To make progress we need to identify the degrees of some of 
the modes either from multi-colour photometry or spectroscopy. 

Key words: stars: individual: EX Cnc, EW Cnc, stars: blue stragglers, stars: variables: 
S Scuti, open clusters: individual: M67 (NGC 2682) 



1 INTRODUCTION 

We present the asteroseismic analysis of a new data set compris- 
ing photometry of the old open cluster M67 (NGC 2682) from nine 
telescopes during a 43-day campaign with a total of 100 telescope 
nights (Stello et al. 2006; hereafter Paper I). The main goal of the 
campaign was to detect oscillations in the stars on the giant branch 
(Stello et al. 2007; hereafter Paper II). Here we analyse the vari- 
ability of the blue straggler (BS) population. 

BS stars are defined as being bluer and more luminous than 
the tum-off stars in their parent cluster, and they are found in all 



open and globular clusters where dedicated searches have been 
made. BSs are important objects since they are directly linked to 
the interaction between stellar evolution of binaries and the clus- 
ter dynamics. In the cores of globular clusters the stellar density is 
high and direct s tellar collisions may e xplain the existence of some 
of the BS stars toavies&Benzlll995l) . However, in open clusters 
the density of stars is muc h lower and stellar collisions are rare 
dMardling & Aarsethl200lh . BSs in open clusters are therefore gen- 
erally thought to be formed as the gradual coalescence of binary 
stars feane t al. 2006). iHurlev et al] J2005h made a simulation of 
M67 that took into account stellar and binary star evolution, in ad- 
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Figure 1. The colour-magnitude diagram of M67 based on standard B — V 
and V magnitudes. The solid line is an isochrone for an age of 4 Gyr. The 
locations of EW Cnc and EX Cnc have been marked with boxes and six 
other blue stragglers in the instability strip (dashed lines) are marked with 
filled circles. 



dition to the dynamical evolution of the cluster. Their simulation 
showed that around half of the BSs were formed from primordial 
binary systems and the other half from binary stars that had been 
perturbed by close encounters with other stars. 

While several of the known BSs in clusters are located near 
the main sequence inside the c lassical instability str ip, only around 
30% show 8 Scuti pulsations ( iGilliland etaDl 19981) . These objects 
are particularly interesting because their masses can in principle 
be inf erred from comparison with pulsation models. iBruntt et al.l 
d200ll) studied the variable BSs in the globular cluster 47 Tuc, and 
could determine the mass of one of them to within 3.5% from a 
comparison with models in the Petersen diagram. 

We have made a search for oscillations in eight BS stars in the 
instability strip in M67 (Fig. [TJ. In particular, we have studied the 
two known variable BS stars EW Cnc and EX Cnc. The 8 Scut i 
pulsations in these stars were discovered by IGilliland et alj J 1 99 lh 
and pu lsation modelling was done in detail by IGillilan d & Brown 
(1992). The current data set has a sign ificantly higher signal-to- 
noise (S/N) than the previous studies (Gilliland & Brownlll992l ; 
Zhang et al. 2005) and a superior spectral window. We have for the 
first time unambiguously detected a very large number of frequen- 
cies, comparable in number to the r ecent ambitious cam paigns on 
field 8 Scuti stars like FG Virginis tereger et alj|2005h . We char- 
acteri ze EW Cnc and EX Cnc u sing a grid of seismic models, but 
unlike lGilliland & Brown! Jl992h we also include the effects of ro- 
tation (in both equilibrium models and seismic oscillations). 

Rotation will affect the internal structure of the stars, in- 
duce mixing processes jZ ahn 1992; iHeger. La nger & Wooslevl 
2000; iReese. Lignieres & Ri eutord 2006), and cause significant 
asymmetric splitting of multipl e ts even in these moderately 
rapid rotating stars dSaiol 1 19811; iDziembowski & Goodd [l992l; 
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Figure 2. Finding chart for the blue stragglers we have analysed with ID 
numbers from Sanders (1977). The image is from the STScI Digitized Sky 
Survey. 



2 BLUE STRAGGLER TARGETS 

The locations of EW Cnc and EX Cnc in the colour-magnitude 
diagram of M67 are given in Fig. Q] We also analysed six other 
BS stars inside the instability strip and they are marked by filled 
circles. The B — V colours and V m agnitudes are taken from 
iMontgomerv. Marschall & Janesl dl993l) . Only stars observed dur- 
ing the present campaign are plot ted. We have plotted a n isochrone 
taken from the BaSTI database jPietrinferni et alj|2004l) for an age 
of 4.0 Gyr with composition Z = 0.0198, Y = 0.2734, adopt- 
ing a distance modulus of (m — M)y = 9.7 . The instability strip 
is mar ked by dashed lines and is taken from iRodrfguez & Breged 
j200ll) . It has been transforme d from the S trbmgren to the Johnson 
system using the calibration in lCoxl (2000). A finding chart show- 
ing the locations the BS targets is given in Fig. [2] 



3 OBSERVATIONS 

The data were collected in 2004 from January 6 to February 17 
during a multisite observing campaign with nine 0.6-m to 2.1-m 
class telescopes. The photo metric reduction was done using the 
MOMF photometric package dKieldsen" &Frandsenlll992l Sects. 5- 
6). The details of the observations and photometric data reduction 
were given in Paper I. 

The choice of filters was made to optimize the signal-to-noise 
for the observations of the giant stars in M67 (see Paper II). All sites 
used the Johnson V filters except Kitt Peak, where a Johnson B fil- 
ter was used. The bright BS star EX Cnc was saturated in about half 
of the B images. We have a total of 450 hours of observation and 
about 15 000 data points for each star. The typical point-to-point 
scatter for EW Cnc and EX Cnc varies greatly f rom site to site, and 
lies in the range 1-8 mmag. By comparison, Gilliland & Brown! 
( 1992) used 37 hours of observation to obtain 2458 CCD frames 
with 2 mmag point-to-point noise, which is similar to our best sites. 

The complete light curve of EW Cnc is shown in the top panel 
in Fig. [3] The middle and bottom panels show the details of the 
light curves of EW Cnc and EX Cnc during 24 hours. We have used 
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Figure 3. The complete light curve of EW Cnc is shown in the top panel. The details of 24 hours with nearly full coverage are shown in the two bottom panels 
for EW Cnc and EX Cnc where different symbols are used for each observing site. The continuous curves are the fits to the light curves. 
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Table 1. For case A and B we give the noise level in the amplitude spectrum 
at high frequencies (760 ± 60 /iHz) and the height of the first alias peak in 
the spectral window. The results for our preferred cases are printed in italics. 





Case A 


CaseB 




^760 


Alias 


A760 Alias 


Star 


|>mag] 


[%] 


[Atmag] [%] 


EWCnc 


55 


23 


37 48 


EXCnc 


48 


21 


36 47 



different symbols for the five sites contributing to the light curve, 
as indicated in Fig. [3] The abbreviations for the five observatories 
are the same as in Paper I. The solid curves are the final fits to the 
complete data set (see Sect. l4.4t . 



4 TIME-SERIES ANALYSIS 

4.1 Data point weights 

As can be seen qualitatively in Fig. [3] and as described in detail in 
Paper I the quality of the photometry varied a lot from site to site, 
in many cases from night to night, and sometimes also during the 
night at a given site. To optimally extract the pulsation frequencies 
from the time series we have assigned weights to each data point. 

We calculated two kinds of weights, "scatter weights" (Wscat) 
that are based on the scatter in the time series and "outlier weights" 
(Wout), which suppress extreme data points. The weights were cal- 
culated for each telescope and each star. The scatter weights were 
determined from the average point-to-point scatter in a large en- 
semble of reference stars, as described in detail in Appendix IA1I 
To calculate the outlier weights, we used the deviation from zero 
of both the target star and the reference stars, as described in Ap- 
pendix |A2] 

4.2 Combining observations in B and V 

At Kitt Peak we used a Johnson B filter while all other sites used 
the Johnson V filter. The Kitt Peak data w ere collected simul tane- 
ously with several other sites (see Fig. 2 in lStello et alj d2006t) ) but 
since the precision is far superior we decided to include it. The B 
filter data comprise 1459 data points (9%) for EW Cnc and 598 
points (4%) for EX Cnc. 

When combining the data we must correct for the dependence 
of the pulsation amplitude on the filter by scaling the B filter data. 
The scaling was found by measuring the amplitude ratios of the 
seven dominant frequencies in EW Cnc when using the B and V 
filters alone. The mean amplitude ratio was 0.74 ± 0.o£J which 
is the scaling we applied. We also scaled the scatter weights found 
in Appendix lAll to ensure that the signal-to-noise is preserved. The 
mean phase shift is 5.7 ± 4.3°, which is not significantly different 
from zero. 

To make sure that our results are not affected by the inclu- 
sion of the scaled B data, we repeated the analysis described below 

1 iKoen et alj 1200 ll) observed two single-mode <5 Scuti stars (HD 199434 
and HD 21190) in the B and V filters and found amplitude ratios of 
0.70 ± 0.01 and 0.70 ± 0.03, respectively. We note that the empirical 
scaling is close to the ratio of the central wavelengths, ie. \g /Xy = 
434 nm/ 545 nm~ 0.80. 



Table 2. Frequency, amplitude, phase, and S/N of frequencies detected in 
EW Cnc for case B. 



ID 


/ I/^Hz] 


a [mmag] 


4> 


S/N 


/1 


217.993(4) 


2.66(2) 


0.05(1) 


48.0 


h 


233.789(6) 


1.79(3) 


0.95(2) 


31.8 


h 


363.811(4) 


1.29(3) 


0.00(1) 


26.9 


h 


322.244(5) 


1.17(3) 


0.76(1) 


23.8 


h 


435.39(1) 


0.94(4) 


0.19(3) 


20.8 


h 


254.66(1) 


0.93(3) 


0.94(3) 


16.2 


h 


308.81(1) 


0.75(2) 


0.66(3) 


14.7 


h 


435.82(3) 


0.55(4) 


0.62(7) 


12.2 


h 


367.83(1) 


0.58(3) 


0.93(3) 


12.2 


ho 


508.98(2) 


0.53(2) 


0.15(4) 


11.6 


hi 


443.07(3) 


0.46(3) 


0.78(7) 


10.1 


hi 


299.25(1) 


0.51(2) 


0.96(4) 


9.7 


ha 


470.16(2) 


0.41(3) 


0.93(5) 


9.1 


fa 


474.36(2) 


0.40(3) 


0.34(5) 


8.8 


/is 


513.75(2) 


0.40(2) 


0.50(5) 


8.6 


/if, 


184.69(3) 


0.46(3) 


0.47(7) 


8.6 


fn 


485.12(2) 


0.38(3) 


0.76(6) 


8.5 


/is 


442.59(3) 


0.37(4) 


0.64(7) 


8.2 


/l9 


530.87(3) 


0.36(2) 


0.62(7) 


7.5 


/20 


375.96(3) 


0.33(3) 


0.09(7) 


7.0 


/21 


354.22(3) 


0.33(3) 


0.92(9) 


6.9 


f'22 


351.48(7) 


0.33(4) 


0.00(6) 


6.8 


j'2'A 


202.14(3) 


0.37(3) 


0.06(7) 


6.7 


f'24 


441.59(2) 


0.30(3) 


0.98(6) 


6.7 


f'23 


557.87(3) 


0.32(3) 


0.15(8) 


6.7 


/2(i 


445.43(3) 


0.29(2) 


0.23(9) 


6.5 


127 


593.32(3) 


0.30(2) 


0.82(7) 


6.2 


f28 


506.27(4) 


0.28(2) 


0.72(9) 


6.0 


f29 


311.87(3) 


0.29(2) 


0.13(8) 


5.7 


ho 


590.05(4) 


0.27(2) 


0.5(1) 


5.7 


hi 


586.55(5) 


0.26(3) 


0.5(1) 


5.5 


hi 


337.43(4) 


0.25(3) 


0.8(1) 


5.2 


h-i 


329.42(2) 


0.25(3) 


0.95(7) 


5.1 


/34 


586.15(5) 


0.24(3) 


0.1(1) 


5.1 


/35 


400.09(5) 


0.24(3) 


0.4(1) 


5.1 


/36 


235.16(4) 


0.28(3) 


0.3(1) 


4.9 


f:w 


475.67(5) 


0.20(2) 


0.5(1) 


4.5 


/38 


413.53(6) 


0.20(3) 


0.7(1) 


4.4 


h'.t 


625.85(4) 


0.21(2) 


0.8(1) 


4.3 


ho 


315.65(4) 


0.20(2) 


0.3(1) 


4.1 


hi 


421.70(4) 


0.19(2) 


0.8(1) 


4.1 



when using only the V data. For both stars we found exactly the 
same frequencies except for two frequencies in EW Cnc and three 
frequencies in EX Cnc. The most significant of these frequencies 
has S/N= 5.2 while the others have S/N below 4.4. The slight 
differences we find in the frequencies and amplitudes are within 
the uncertainties found in Appendix [B] We are therefore confident 
about combining the V data with the scaled B data. 

4.3 The optimal amplitude spectrum 

We calculated the amplitude spectrum both with the optimal spec- 
tral window and with the lowest noise level (see Kjeldsen et al. 
2005). To do this we made the time-series analysis for two differ- 
ent approaches: 

Case A: We computed the average light curve by binning all data 
collected within 8-minute intervals. In the binning process we took 
into account the different quality of the data points by using weights 
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Figure 4. Amplitude spectra of EW Cnc (left panels) and EX Cnc (right panels) for the binned light curves (case A, top panels) and the noise optimized light 
curves (case B, bottom panels). The residual spectrum is shown below each plot and the insets show the spectral windows. In the residual spectra the dashed 
grey lines mark the average noise level. 



Wscat • Wont for n — 2 in Eq. IA3I We did not use weights when 
computing the amplitude spectrum a second time. 
CaseB: No binning, but each data point was given the weight 
Wscat ■ Wont when calculating the amplitude spectrum. 



It was case A that provided the optimal spectral window. The 
amplitude spectra of EW Cnc and EX Cnc for case A are shown in 
the top panels in Fig. [4] The residual spectrum after extracting the 
oscillation frequencies (see Sect. |4.4t is shown in each panel and 
the insets show the spectral windows, which have sidelobes at ±1 
c d" 1 and ±2 c d" 1 of ~ 23%. 

The lowest noise level was found in case B. We tried both 
n = 1 and 2 in Eq. |A3| but found that for n = 1 the noise level 
was 5-10% lower, and in addition the ±1 cd _1 sidelobes were 
slightly lower. The fact that we found the lowest noise level for 
n = 1 indicates that l/(H4cat ■ Wout) is the best approximation 
to the true variance in the time series. We note that this was also 



concluded bv lHandlei] J2003h . The amplitude spectra for case B are 
shown in the bottom panels in Fig. [4] 

While the noise level is significantly lower for case B, the 
spectral window has much lower sidelobes in case A. In Table Q] 
we compare the noise levels and aliases. The noise was calculated 
in a frequency band of 760 ± 60 /iHz in the residual spectra. 

In Appendix |B] we have made simulations of the time series 
to check whether case A or B gives the most robust results. The 
conclusion is that case A is preferred for EX Cnc due to the high 
number of frequencies found in a relatively narrow interval, while 
we use case B for EW Cnc. The results for these cases are printed 
in italics in Table [T] 

4.4 Fourier analysis of EW Cnc and EX Cnc 

For each site the light curves were high-pass filtered to remove slow 
trends. In this way we suppressed variations with frequencies be- 
low 80 /iHz, corresponding to periods longer than ~ 3.5 hours. 
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Table 3. Frequency, amplitude, phase, and S/N of frequencies detected in 
EX Cnc for case A. 



ID 


f ImHz] 


a [mmag] 


4> 


S/N 


fi 


226.910(4) 


3.87(6) 


0.94(1) 


36.2 


h 


238.883(5) 


3.74(6) 


0.60(1) 


35.0 


h 


240.297(8) 


2.31(7) 


0.97(2) 


21.7 


h 


191.464(9) 


2.08(7) 


0.45(2) 


18.3 


h 


226.45(1) 


1.79(8) 


0.28(3) 


16.7 


h 


205.49(1) 


1.68(8) 


0.14(3) 


15.2 


h 


228.66(2) 


1.33(8) 


0.57(4) 


12.4 


h 


215.33(1) 


1.34(8) 


0.41(4) 


12.3 


h 


196.96(1) 


1.18(7) 


0.71(4) 


10.5 


iio 


190.92(1) 


1.16(8) 


0.46(3) 


10.2 


in 


217.43(2) 


1.05(8) 


0.28(5) 


9.6 


/l2 


211.82(2) 


1.04(7) 


0.77(4) 


9.5 


/l3 


196.15(2) 


0.83(7) 


0.03(4) 


7.4 


ii4 


193.93(2) 


0.79(7) 


0.11(6) 


7.0 


fir, 


219.95(2) 


0.75(7) 


0.89(6) 


6.9 


fid 


236.34(2) 


0.60(5) 


0.04(7) 


5.6 


fl7 


203.19(3) 


0.62(7) 


0.64(8) 


5.6 


flS 


223.01(4) 


0.59(7) 


0.19(9) 


5.5 


fw 


258.86(3) 


0.58(6) 


0.27(8) 


5.5 


ho 


199.45(3) 


0.59(7) 


0.76(7) 


5.3 


hi 


234.39(4) 


0.55(6) 


0.5(1) 


5.2 


f-22 


229.14(5) 


0.55(5) 


0.5(1) 


5.2 


,hi 


250.65(2) 


0.48(8) 


0.85(6) 


4.5 


hi 


256.86(3) 


0.47(6) 


0.55(9) 


4.4 


hr, 


182.09(3) 


0.50(6) 


0.69(8) 


4.4 


he 


148.54(3) 


0.49(6) 


0.20(7) 


4.3 



The pulsation frequencies are found mainly in the intervals from 
200-600 /iHz and 150-350 /iHz for EW Cnc and EX Cnc resp ec- 
tively. We used the PERIOD04 package bv lLenz & Bregen ( 120051) to 
fit the light curves in the frequency range 100-1000 /iHz, as fol- 
lows. The amplitude spectrum was calculated using the light curve 
from case B for EW Cnc and case A for EX Cnc (see Sect. |4.3I >. 
The highest peak was then selected and the frequency, amplitude 
and phase were fitted. This was done iteratively while in each step 
always improving the frequencies, amplitudes and phases of previ- 
ously extracted frequencies. 

We extracted 41 and 26 frequencies with S/N > 4 in the two 
stars, and the parameters are given in Tables [2] and [3] The uncer- 
tainties are based on simulations, as described in Appendix |B1 To 
determine the S/N of the extracted frequencies, we estimated the 
noise level in the amplitude spectrum around each frequency. This 
was done after cleaning the amplitude spectra so only peaks with 
S/N < 3 remained. The noise was calculated in bins with a width 
of 35 /iHz and the noise estimate is the average of the three near- 
est bins for each frequency. The noise levels are marked by grey 
dashed lines in the residual spectra in F ig. [4] 

In comparison iGilliland & BrowrJd 19921) detected 10 frequen- 
cies in EW Cnc and 6 frequencies in EX Cnc. We recovered these 
frequencies except one of their wea ker modes at 423. 2 /iHz in 
EW Cnc and at 186.5 /iHz in EX Cnc. lZhang et alj d2005t) detected 
4 frequencies in EW Cnc and 5 frequencies in EX Cnc. They used 
data from a single observatory, and this explains why some of their 
frequencies are offset from ours by ±1 c d _1 . 

To check whether the frequencies we detected in the time- 
series analysis were reliable, we made the time-series analysis for 
both cases A and B. In both EW Cnc and EX Cnc we found that 



nearly all frequencies with S/N > 7 were detected in both cases, 
while this was not true for the weaker frequencies. Further, for 
the weaker frequencies that were detected in both cases there was 
sometimes an offset of 1 c d~ .In those cases, we chose the solu- 
tion for case A, which has the better spectral window. As an addi- 
tional check we used the 20 days in the time series with the best 
coverage and were able to recover all frequencies with S/N > 6. 
We also used only the data set from La Silla (14 consecutive 
nights of observation) but found that for the weaker frequencies, 
the ±1 c d _1 alias was often identified as the most significant fre- 
quency, which is not surprising given the high sidelobes from the 
single-site time series. 

Several of the 41 frequencies found in EW Cnc are found to 
be linear combinations; examples are /s ~ /13 ~ 2/i and fig ~ 
fs — f'2 — fi2 — hs- We searched for combinations within twice the 
resolution of the data set which is A/ = 2/T b s = 0.63 /iHz. The 
condition is expressed as fi ~ nfj + mfk for any integer value 
n,m £ [—2; 2] and i,j, k £ [1; 41] and we found 13 combina- 
tion frequencies. We compared this with 1 000 sets of 41 randomly 
distributed numbers in the same interval as the observed frequen- 
cies, and we found 18 ± 6 combinations. Hence our detection of 13 
linear combinations in EW Cnc is not surprising. However, none 
of the frequencies found in EX Cnc are found to be linear com- 
binations. This is also the case for 1 000 simulations of randomly 
distributed frequencies. The reason for the different results for the 
two stars is that the frequencies in EX Cnc are found in a much 
narrower frequency interval than for EW Cnc. 



4.5 Interpretation of the EX Cnc spectrum 

In both cases A and B, we see a clear pattern of nearly equally 
spaced peaks for EX Cnc (FigQ. In the top panel in Fig.f5]we show 
the amplitude spectrum of EX Cnc (case A) in more detail. To guide 
the eye we indicate a comb pattern with a separation of 1 cd~* 
(11.574 /iHz) with dashed vertical lines. In the bottom panel the 
frequencies detected in case A (listed in Table [3} are marked by 
vertical lines. 

We note that each dominant peak in the comb pattern con- 
sist of several closely spaced peaks with remarkable similarity to 
what is seen for stochastically excited and damped oscillations, 
also called solar-like oscillations, w hich are driven by convection 
dAnderson, Duvall & Jefferieslll990l) . We investigated whether the 
observed regular pattern and the closely spaced lower amplitude 
peaks around each dominant peak could be explained by 

(i) Stochastically excited and damped oscillations of equally 
spaced modes. 

(ii) The spectral window (alias peaks). 

(iii) Closely spaced modes. 

(i) In 5 Scuti stars, the oscillations are caused by the n- 
mechanism (opacity driven) and are believed to be coherent, 
which will produce a single isolated peak for each oscillation 
mode in the amplitude spectrum. Although damped oscillations 
driven by convection might be present in some S Scuti stars 
dSamadi. Goupil & Houdeki |2002|) . the expected frequencies and 
amplitudes are very different from what we see in Fig. [5] which 
is consistent with opacity driven oscillations. The calculations by 
ISamadi. Goupil & Houdeki d2002h show frequencies in the range 



2 The S/N was measured for case B. 
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Figure 5. The top panel shows the amplitude spectrum of EX Cnc for case A. The inset shows the spectral window on the same horizontal scale as the main 
Figure. The vertical dashed lines are separated by 1 c d~ 1 (11 .574 pHz). In the bottom panel the detected frequencies are marked by vertical lines. 



Table 4. Properties of the light curves and amplitude spectra of the observed BS stars in M67. Star ID from Sanders (1977) and ID, V, and B — V from 
Montgomery, Marschall, & Janes (1993) are given. N is the number of data points. The point-to-point scatter at the LaS and LOAO sites and the noise level in 
the amplitude spectrum in the range 350 ± 150 pHz and 760 ± 60 pMz are listed. The last column contains membership probabilities from Sanders (1977). 















°"LaS 


"LOAO 


^350 


^760 


■^Member 


Name 


Sanders 


MMJ 


V 


B — V 


N 


[mmag] 


[mmag] 


[>mag] 


[/imag] 


[%] 




S1466 


6511 


10.60 


0.34 


13316 


1.7 


3.3 


66 


41 







S1434 


6510 


10.70 


0.11 


3312 




3.3 


372 


160 


91 


EX Cnc 


S1284 


6504 


10.94 


0.22 


14660 


1.9 


2.9 


240 


65 


95 




SI 066 


6490 


10.99 


0.11 


16954 


1.7 


2.6 


46 


35 


90 




SI 263 


6501 


11.06 


0.19 


17408 


4.6 


2.7 


89 


69 


89 




S968 


6479 


11.28 


0.13 


9851 


1.5 


2.9 


55 


37 


95 




S752 


6476 


11.32 


0.29 


3316 




2.8 


194 


120 


95 


EWCnc 


S1280 


5940 


12.26 


0.26 


15524 


1.9 


3.3 


188 


46 


93 



400-1200 /iHz and amplitudes of roughly 100 ppm for their se- 
lected models. These values are consistent with what we found 
specifically for this star using the parameters in Table[5] a mass of 
2 M©, and the usual scalin g relations to predict the characteristics 
of so lar-like oscillations dKieldsen & Bedding! 1 19951 : iBrown et alj 
Il99ll) . We find the central frequency (maximum amplitude) to be 
~ 630 fiHz, amplitudes 8L/L ~ 40 ppm, and mean fre- 
quency separation Auo — 40 /xHz). These properties of solar-like 
oscillations are clearly inconsistent with the observed amplitude 
spectrum. 

(ii) Simulations showed that the regular pattern was unlikely to 
be the result of only two or three modes. To reproduce the observed 
amplitude spectrum required at least 4-5 modes, equally spaced 



by roughly led -1 . The low alias sidelobes, especially in case A, 
imply that the pattern is intrinsic to the star. Many of the lower 
amplitude peaks around each dominant peak seem to be separated 
by exactly 1 or 2cd _1 from a neighbouring peak. However, the 
strength of these lower amplitude peaks are in many cases higher 
than expected from the strength of the sidelobes in the spectral win- 
dow. 

(iii) Closely spaced peaks are not uncommon in 5 Scuti stars 
( iBreger & Bischo f 2002). B reger & Pamvatnvkhl j2006h noted that 
18 pairs (A/ < 1 /xHz) exist in the S Scuti star FG Virginis 
jBreger et al.ll2005T) . They found evidence for the frequency dou- 
blets being real. As seen in Fig. [5] we also detect several close 
pairs. 
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In summary, damped oscillations driven by convection seems 
to be a very unlikely cause for the characteristic pattern we see in 
Fig. [5] We believe the pattern is most likely due to a nearly reg- 
ular series of at least five modes (excited by the k mechanism) 
which, due to a general spacing very close to 1 c d _1 , shows a large 
number of small peaks around each dominant peak caused by the 
spectral window. This is further enhanced by the presence of lower 
amplitude modes close to the dominant mode s (see bottom panel 
in Fig . Our interpretation is supported by iGilliland & Brownl 
dl992h who found frequency separations of ~ 5 /iHz and ~ 10 /iHz 
in this star, and suggested that the latter might be due to rotational 
splitting. 



4.6 A search for pulsations in six other BS stars 

In addition to EW Cnc and EX Cnc, we observed six other BS stars 
in the instability strip as shown in Fig.[T] To search for pulsations 
in the stars we calculated the amplitude spectra in the range 0- 
1000 /.iHz, using weights as for case B (Sect. [4~3b . In Table [4] we 
summarize the properties of the BS stars observed during the cam- 
paign. For completeness, EW Cnc (S 1280) and EX Cnc (SI 284) are 
incl uded in the Table. We give the ID numbe rs from lSandersI J 19771) 
and lMontgomerv. Marschall & Janes] dl993l) . and the V magnitude, 
and B — V colour from the latter. The total number of data points 
N and o int (Eq.|AT) for the La Silla and LOAO (Mt. Lemmon Op- 
tical Astronomy Observatory) sites indicate the quality of the data. 
We note that La Silla had very good weather conditions during the 
campaign, while LOAO had worse conditions (for more details see 
Paper I). The mean level in the amplitude spectrum in the range 
where oscillations are expected is denoted A350 (200-500 ^Hz), 
while that at slightly higher frequencies j4.76u (700-820 ^iHz) is a 
measure of the white no i se. In the last column, the membership 
probability from Sanders 1 1977b is given. W e no te that more recen t 
radial velocity studies bv lGirard et all d 19891) and lZhao etal1dl993h 
are in general agreement with this. 

Two of the BS stars (S1434 and S752) were only observed by 
a few sites. Hence, for these stars the noise level in the amplitude 
spectrum is about a factor 4-6 higher than for the other stars. We 
find no significant peaks in the range 0-1000 /iHz above S/N of 4f| 

in any of the BS stars except EW Cnc and E X Cnc. 

Our results are in agreement with Sandquist & Shetronel 

J2003I) who looked for photometric variability in several BS stars 
in M67. Except for EW Cnc and EX Cnc they found no S Scuti 
pulsations in their sample, which includ ed S752, S968, S1066, and 
SI 263. The latter was also observed by IGilliland & Brownl fl992T) 
who placed an upper limit of variation at 0.2 mmag (correspond- 
ing to ~ 5 mmag scatter in the time series) for frequencies above 
130«Hz, in agreem ent with this s t udy a nd Sandquist & Shetronel 
(2003). We note that lStassun et al.l d2002l) reported unusually high 
scatter of up to 30 mmag in this star from a photometric B VI time- 
series study. The relatively high noise level for SI 263 compared to 
the other BSs in our data (see OLaS in Table [4]l is due to the star 
being blended (see Fig.[2]l. 



Table 5. Fundamental parameters of EW Cnc and EX Cnc. 



5 THEORETICAL MODELS OF EW & EX CNC 

To mod el EW Cnc an d EX Cnc we used the stellar evolution code 
CESAM dMorell 19971) . The effects of rotation are taken into account 



3 Defined as 4 • A350 where A350 is from Table|4] 





EWCnc 


EX Cnc 




T efi IK] 


log 9 


T eff [K] log g 


Mathys et al. 1991 


8090 


4.33 


7750 3.79 


Gilliland & Brown 1992 


7960 


4.20 


7900 3.88 


Landsman et al. 1998 


8090 


4.12 


7610 3.78 


uvby (this study) 


7800 




7840 


Geneva {this study) 


8080 


4.59 





Table 6. Parameters of theoretical pulsation models for EW Cnc and 
EX Cnc, which are inside the photometric error boxes shown in Fig. \§\ 
For each mass the range of effective temperatures is given. For each star we 
computed models with two different rotational velocities given in the last 
column. 



Star 


M/Mq 


T eff [K] 


Vrot [kms 1 ] 


EWCnc 


1.7 


7930-8110 


90, 150 




1.8 


7930-8210 






1.9 


7930-8210 






2.0 


7940-8200 




EX Cnc 


1.8 


7450-7620 


65, 120 




1.9 


7450-7740 






2.0 


7480-7740 






2.1 


7470-7710 






2.2 


7480-7730 






2.3 


7630-7740 





in the equilibrium equations as a first-order perturbation to the local 
gravity, while assu ming rotation as a rigid body. The CEFF eq uation 
of state was used dChristensen- Dals gaard & Daepperi Jl992D and 
OPAL opacities were adopted Jlglesias & Rogers! Il996t) . A more 
detailed description of the input physics for the CESAM code was 

given bv lCasasetaT] ( l2006l) . 

We used the oscillation code FILOU dTran Minh & Leonl 19951: 
ISuarez!l2002h to calculate the oscillation frequencies from the evo- 
lution models. We assume rigid-body rotation and the effects up 
to second order in th e centrifugal and Coriolis forces w ere in- 
cluded. We note that iReese. Lignieres & Rieutordl |2006) found 
that this perturbation approach may be invalid even at moderate 
rotational velocities. Effects of near degeneracy are expected to 
be significant (Soufi et al. 19 98) as was described in detail by 
ISuarez. Goupil & Morell d2006l) . 



5.1 Fundamental atmospheric parameters 

The fundamental a tmospheric para mete rs of EW Cnc and EX Cnc 
were estimated by iMathyj dl99ll) a nd IGilliland & Brownl J 19921) 
based on Stromgren photometry by iNissen. Twarog. & CrawforT 
(1987) but applying differ ent calibrations dMoon & Dworetsk 
Il985l: iPhilip & Reive j[l97l). We used the more recent calibration 
by iNapiwotzki. S choenber ner & Wenskel dl993h . while assuming 
the mean cluster reddening from Nissen et al. (198 7). For EW Cnc 
we used Gene va photometry from Rufener ( 1976 ) an d appli ed the 
calibration by iKunzli et alj dl9971) . lLandsman et alj d 19980 con- 
strained Teg of both stars based on UV measurements. We summa- 
rize the results in Table [5] Typical uncertainties on T c g and log g 
are 150 K and 0.2 dex and there is good agreement between the 
different calibrations. 
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Figure 6. The log g— T e ff diagram showing selected evolution tracks from 
our model grid. The 1 — a photometric error boxes of EW Cnc and EX Cnc 
are indicated as grey squares. 



An important limitation of our modelling is the uncertainty 
on the composition of the stars. The composition will depend on 
the formation scenario of these BS stars, whic h we certainly can- 
not be gin to constrain with the current data set. lGilliland & Brownl 
d 19921) discussed different formation scenarios for the same stars 
from either a direct collision or gradual coalescence. They made 
the important point that it would be difficult to detect the ef- 
fects of a 5% increase of helium, eg. as a result of mixing two 
evolved main sequence stars. For simplicity we assume that the 
metallicity of EW Cnc and EX Cnc is t he same as for the cluster. 
iMontgomerv. Marschall & Janes! i 19931) found [Fe/H] = -0.05 
but they did not estimate the uncertainty; a conservative estimate 
is ±0.10 dex. 



5.2 The model grid 

We have calculated a grid of stellar models consistent with the 
fundamental parameters of EW Cnc and EX Cnc. A few evolu- 
tion tracks are shown in Fig. [6] in the loggr versus T c a -diagram, 
where the solid and dashed lines correspond to rotational veloci- 
ties of 65 and 120km s -1 . The photometric error boxes are shown 
as grey squares. The ce nter of each box is the T e s and log g from 
lLandsman et alj Jl998h and the half-widths are equal to the esti- 
mated uncertainties of 150 K and 0.2 dex. We have computed pulsa- 
tion frequencies for the models inside the photometric error boxes. 
The mass ranges are 1.7-2.0 Mq for EW Cnc and 1.8-2.3 Mq for 
EX Cnc in steps of 0.1 Mq. Further details of the grid of pulsation 

models are given in Table|6] 

Analyses of stellar spectr a (Peterson, Carnev & Latham 1984; 
Manteiga. Martinez Roger & Pickles! 1 19891 ; IPritchet & Glaspevl 



19911) reveal moderate projected rotational velocities for EW Cnc 



and EX Cnc, and we have used average values of v sin i — 90 
and 65kms~ 1 , respectively. However, the inclination angle is 
unknown, and therefore we calculated models and corresponding 
oscillation frequencies for the minimum and a moderate value 
of v TOt for each star. The values in the last column in Table [6] 
correspond to inclination angles of i — 90° and 45°. 



5.3 Comparison of observed and theoretical frequencies 

In Fig. [7] we show the number of modes as a function of frequency, 
in bins of 50 /iHz. The black histograms are for the observed fre- 
quencies for EW Cnc (left panels) and EX Cnc (right panels), 
which are listed in Tablesf2]and[3] The solid histograms are for fre- 
quencies from pulsation models and the colours indicate how high 
angular degree, I, we have included^. The top and bottom panels 
are for models with high and low mass. The low-mass models both 
have a mass 1 . 8 Mq , while the high-mass models correspond to the 
evolution tracks that graze the photometric error boxes, ie. 2.0 and 
2.3 Mq for EW Cnc and EX Cnc, respectively. We used the T e ff 
adopted for the centres of the photometric error boxes in Fig. [6] 

By inspecting the histograms for the models in Fig.|7]one can 
see that the number modes increase when increasing the mass or 
decreasing T e g . Histograms for models with different rotational ve- 
locities are nearly identical since the frequency shifts are relatively 
small (~ few /iHz) for these moderate velocities. The histograms 
of all pulsation models within the photometric error boxes are in 
agreement with the observations, but only when including modes 
with degrees I = 0, 1, and 2. However, if the masses are lower than 
1.7 and 1.8 Mq for EW Cnc and EX Cnc, respectively, it will re- 
quire lower temperatures, which are outside the photometric error 
boxes. In other words, if a more accurate T e g becomes available 
we can put a lower limit on the mass. 

In EW Cnc the detected frequencies are evenly distributed 
over a wide interval, while in EX Cnc they are concentrated in 
a narrow band. The horizontal bars labeled "A" in Fig. [7] indicate 
where most observed frequencies are found. In these frequency in- 
tervals we estimate that we have detected 40-60% of the modes in 
EW Cnc and 70-100% in EX Cnc. These results rely on the as- 
sumptions that the considered ranges in mass and T e « are realistic 
and that only I ^ 2 modes have detectable amplitudes. 

To further investigate the distribution of peaks and search for 
recurring patterns we computed autocorrelations for the frequen- 
cies from the observations and the pulsation models. This method 
can in principle be used to extract important quantities for astero- 
seismic modelling like the large and small separation and rotational 
splitting. This would enable us to constrain the mass and temper- 
ature of the stars. Since we have no information about the relative 
amplitudes of the theoretical frequencies, each frequency is repre- 
sented by a Gaussian of unit height and with a FWHM equal to 
twice the formal frequency resolution of the data set. To be able to 
compare the autocorrelations of the observed and theoretical fre- 
quencies, we follow the same procedure for the observations, ie. 
only using the frequencies and not the observed amplitudes. We 
use the frequencies from the regions where the majority of the fre- 
quencies are observed. These regions are marked by horizontal bars 
labeled 'A" in Fig. [7] We tried different widths of these bands, but 
the conclusion summarized here is the same. 

In Fig. [8] we show autocorrelations for theoretical and ob- 
served frequencies for EW Cnc (left panels) EX Cnc (rightpanels). 
The models have the same masses and temperatures as shown in 
the histograms in Fig. [7] In the top panel only radial modes are 
included. The results for EX Cnc shows the large separation at 
~ 23 pHz for the model with high mass, while for the less mas- 
sive model there is a low peak at twice the large separation (see top 



4 In the integrated light from the stars we expect to observe only low degree 
modes. Modes with I ~^ 3 will have very low observed amplitudes due to 
geometrical cancellation effects when averaging over the stellar disc. We 
note that for each mode we include all 21 + 1 possible values of m. 
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Figure 7. Histograms showing the number of frequencies in EW Cnc (left panel) and EX Cnc (rightpanel) in frequency bins of 50 /^Hz. The filled histograms 
are for pulsation models and the colour describes how high angular degree is included. The thick black lines are for the observations. The top and bottom 
panels are for pulsation models with high and low mass with the adopted T c g . The models in the bottom panels both have mass 1.8 Mq, slightly different 
v TO t, and quite different T c ff . The horizontal bar marked "A" marks the region used for the autocorrelation in Sect. 15.31 
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Figure 8. Comparison of the autocorrelations of frequencies from pulsations models and the observed frequencies in EW Cnc (left panels) and EX Cnc 
(rightpanels). In the top panels we included only radial modes while in the middle panel we included modes with I = 0,1 and 2. We show results for the same 
models as in Fig. [7] The bottom panels are the autocorrelations of the observed frequencies. 



right panel in Fig.[8jl. In the middle panel we included all model fre- 
quencies with I = 0, 1, and 2; the higher number of modes makes 
the interpretation of the autocorrelation much more difficult. In the 
bottom panels we show the autocorrelations of the observed fre- 
quencies. There is an intriguing similarity between the autocorre- 
lation of the observed frequencies in EX Cnc and the model with 
high mass (the middle and bottom right panels in Fig.[8jl. However, 



we are extremely cautious not to over-interpret the only marginally 
significant autocorrelation peaks at 2, 6, 12, and 24/iHz. 

Although there appears to be qualitative agreement between 
the observations and the theoretical models, we cannot at this stage 
use the method to determine which combination of mass, evo- 
lutionary stage, or rotational velocity best describes the observa- 
tions. To make progress and make a detailed comparison of the ob- 
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served frequencies with theoretical models would require that we 
have identified at least some of the dominant modes. This can be 
achieved from the study of line-profile variations or from compari- 
son of phase shifts of the frequencies in different photometric filters 
dDaszvriska-Daszkiewicz. Dziembowski & Pamva tnvkh 20031). We 
cannot apply the latter method since most of our observations were 
done in the V band with only a few nights of data in B. 



6 CONCLUSION 

We have analysed a new photometric data set of eight BS stars in 
M67. Only in the two known variable stars EW Cnc and EX Cnc 
do we find S Scuti pulsations. We detect 41 frequencies in EW Cnc 
and 26 frequenci es in EX Cnc. Compared to the previous multi- 
site campaign bv lGilliland & Brown! jl992T) we have detected four 
times as many frequencies in these two stars, and the number of 
frequencies is comparable to the best-studied field S Scuti stars. 

The frequencies with the highest amplitudes in EX Cnc are 
separated by roughly 1 cd _1 (see Fig. [5} and our extensive multi- 
site data are therefore essential to study the pulsations in this star. 
We have made simulations of the light curves which show that the 
frequencies are recovered in all cases (see Appendix [E). We re- 
peated the Fourier analyses on subsets of the data to confirm that 
the same frequencies were extracted. Thus, we claim that the ob- 
served pattern cannot be due to the spectral window, but represents 
a real property of the distribution of frequencies in EX Cnc. Addi- 
tional support for this claim is the fact that a similar pattern is not 
seen in EW Cnc. 

We computed a grid of theoretical pulsation models that take 
the effects o f rotation into account, wh ich were not included in ear- 
lier models ( Gilliland & Brownll 1992b . We find general agreement 
between the observed and computed frequencies when comparing 
their distribution in histograms. We can claim to have detected non- 
radial modes and in the range where most of the frequencies are 
observed, and at least 40% and 70% of all low-degree (/ ^ 2) 
modes are excited in EW Cnc and EX Cnc, respectively. This is 
perhaps the most important result of the current study, since in the 
past decade of observational work on S Scuti stars only a tiny frac- 
tion of the frequencies expected from theoretical models were in 
fact detected. A reasonable explanation is that high S/N and good 
frequency resolutio n is needed to be abl e det ect the weakest mod es 
as was suggested by Breg er et"al] J2005t) and lBruntt etafl J2007I) . 

The autocorrelation of the observed frequencies does not show 
any significant peaks in any of the stars. However, autocorrelations 
of the frequencies from theoretical pulsation models show that if 
indeed modes with degree I = 0,1, and 2 are observed, we cannot 
expect to unambiguously detect the large separation from such an 
analysis. This result was obtained while assuming that all ampli- 
tudes are equal, since no amplitude information is available for the 
theoretical models. 

To improve on the interpretation of the rich frequency spectra 
of EW Cnc and EX Cnc, and other well-observed S Scuti stars in 
general, we see two possibilities: 

• On the theoretical side, it would be very helpful to have esti- 
mates of the amplitudes of the modes. If relative amplitudes could 
be reliably estimated from theory the autocorrelation might be a 
useful method for measuring the separation between radial modes 
and hence constrain the density of S Scuti stars. 

• On the observational side we need to be able to iden- 
tify the degrees of at least the main modes. This can 
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Figure Al. Internal scatter in the light curves of bright stars observed at the 
LOAO site. The grey line marks the fiducial o-i ntj fidu C . The filled circles 
mark the reference stars used for measuring weights. The inset shows the 
location of the stars in the colour-magnitude diagram of M67. 



in principle be done by comparison of phases and ampli- 
tude ratios for different narrow-band filters as has been at- 



tempted for other 8 Scuti stars I Mova. Garrido & Dupretl [20041 : 



Daszvriska-Daszkiewicz et 



st ars m 
ai]l2005h~ 



We have too few observations in the B filter to follow the latter 
suggestion, and multi-filter (eg. Stromgren by) observations should 
be considered for similar campaigns in the future. The points men- 
tioned h ere also apply when planning groun d-based support for the 
CoRo T jMichelet al1l2006t) and Kepler dBasri. Borucki. & Kochl 
2005) satellite missions. These missions will provide photometric 
time series with long temporal coverage of several S Scuti stars, but 
only in a single filter. 

On request we would gladly provide the light curves presented 

here. 
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APPENDIX A: DATA POINT WEIGHTS 

This Appendix contains a detailed description of how we calculated 
weights using an ensemble of reference stars (cf. Sect. 14. it . The 
reference stars were chosen from the most isolated stars, ie. stars 
with a maximum of 10% flux from neighbouring stars within their 
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Figure A2. The top panel shows the light curve of a reference star from 
WFI. The grey points have scatter weights Wscat < 25%. The bottom 
panel shows the scatter weight found using all the reference stars. The 
dashed horizontal line marks the 25% level. 



photometric apertures. Furthermore, we required that their point- 
to-point scatter was close to the theoretical limit as determined by 
photon noise. The observations from LOAO had the greatest field- 
of-view with a total of 358 stars. For this site we used M — 65 
reference stars. For the site with the smallest CCD we used only 
40 reference stars. In Fig. lAll we show the average point-to-point 
scatter (ermt, see Eq. lAll below) for the stars at the LOAO site plot- 
ted versus V magnitude. The solid grey line marks the fiducial 
fflnt.fiduc that represents the typical scatter for the best stars at a 
given magnitude. The reference stars are marked with filled circles, 
and the inset shows the location of the stars in the colour-magnitude 
diagram. 



Al Determination of scatter weights: W BC at 

For each reference star we calculated the internal standard deviation 



JV 

°L = J2( mk ~ ™<.-+i) 2 /[2(iv - 1)], 



(Al) 



where is the magnitude of the fcfh measurement and N is the 
number of data points. For each reference star r and for each data 
point fc, we measured how much the data point itself and the two 
neighbouring data points (k — 1 and k + 1) deviated from zero, ajf] 



|Amk+i| + |Amfc| + Am fc _i| 

3 0"int,fiduc 



(A2) 



where Am is the deviation from zero and <Ti n t,fiduc is the fiducial 
scatter (cf. solid curve Fig.lAfl. Thus, S r is a measure of the actual 
scatter normalized by the expected scatter for a star of that magni- 
tude. 



5 The neighbouring data points were required to be within a certain time 
limit, set arbitrarily to ten times the median time step between frames (for 
a given telescope). If there was only one neighbouring data point, we used 
only this to compute the deviation S r . If there were no neighbouring data 
points within the time limit (rarely the case), we assigned a weight that was 
the mean of the weights of data points just before and after the point in 
question. 
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Figure A3. The light curve of a reference star is plotted versus a target star 
(EW Cnc) for data from WFI. The large black dots have outlier weights 
W r < 0.5. We plot the "3 a ellipse" in grey colour (see text for details). 



Based on all M reference stars, we then computed the average 
deviation of the fcth data point as (5) = y~) 5 r /M, and calcu- 
lated the scatter weight as 



W s < 



(A3) 



«5> + <5floor)" ' 

where we used <5fl oor = 1 mmag to avoid that a few very good 
stars dominate. The optimal choice of the exponent n is discussed 
in Sect. |431 

In Fig. IA2l we show the light curve versus data point number 
for one reference star in the top panel. These data are from the 
Wide-Field Imager (WFI) at Siding Spring, Australia. In the bottom 
panel we plot the derived scatter weights normalized so that the 
highest weight equals 1 . The dashed horizontal line marks the limit 
for data points with scatter weights H4cat < 25%. These points 
are plotted with grey symbols in the top panel as well. 

A2 Determination of outlier weights: W ou t 

We want to assign low weight to data points that deviate signifi- 
cantly from the mean magnitude. We calculated the outlier weight 
for each data point for the rth reference star as 



W r = 





Ami 


b 






+ 








a t ■ a t 




a r ■ o~ r 





(A4) 



and the final outlier weight for the target star is the average value 
W out = Yli^—i W r /M using all M reference stars. In Eq. IA4I 
Arm is the deviation from zero for the target star and at is an esti- 
mated scatter. Likewise, Am r and a r are the corresponding values 
for the rth reference star. The values of a t and a r determine how 
much a data point must deviate before the outlier weight becomes 
significantly smaller than 1. For the target star we used a t = 3 and 
for the reference stars we used a r — 2. Thus, data points that lie 
within 3 a for the target star, and within 2 a for the majority of the 
reference stars will not get low outlier weights. The exponent b de- 
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Figure A4. In the top panel the light curve of a reference star is shown, 
where filled circles have reference star outlier weights W r < 0.5. The light 
curve of the target star (EW Cnc) is shown in the bottom panel, where filled 
boxes have final outlier weights Wout < 0.5 (see text for details). 



scribes how steeply W r decreases from one to zero. We used the 
value 6 = 4. To estimate the scatter a in Eq. lA4l for the target and 
reference stars, we use the scatter weights W sca t from Sect. lAll We 
approximate a — (Tint ■ where <7i n t is computed from Eq. lAll and 
V = (WscaO/Wscat, where (W S cat) is the average scatter weight 
over the entire series. Thus, r\ is a measure of the relative change in 
scatter from frame to frame. 

To illustrate the process we have plotted the light curve of a 
reference star versus the target star EW Cntjf] in Fig. [A3] The grey 
ellipse has half-major and half-semi-major axes corresponding to 
3 (Tint for each star. Intuitively, there are too many points outside 
this "3 a ellipse". If the noise were Gaussian, the fraction of points 
that would fall outside the 3a limit would be only P>3 CT (2D) = 
1.1%. This corresponds to only 24 out of 2150 data points from 
WFI, while about 100 are found in the data. The large black dots 
in Figure [A3] have outlier weights W r < 0.5 and small points have 
1 ^ W r > 0.5. 

In Fig. lA4l we show the light curves of EW Cnc and the same 
reference star as in Fig. IA3I For the reference star (top panel) 
the filled circles mark data points with low outlier weights, ie. 
W r < 0.5. The data points of the target EW Cnc are shown in 
the bottom panel in Fig. lA4l We have marked the points with final 
outlier weights Wout < 0.5 with filled squarefl It is interesting to 
note that several data points are not obvious outliers. In fact some 
of the marked outliers have eg. Am < 2a, where a is the local rms 
scatter. The reason for these points having Wout < 0.5 is that they 
were outliers in most of the reference stars. 



APPENDIX B: LIGHT CURVE SIMULATIONS 

Simulating the light curves has two goals: 



6 Note that the stellar oscillations were subtracted before calculating out- 
lier weights. 

7 Recall that Wout is the mean of W r for all reference stars r. 



• To measure our ability to extract the inserted modes for 
cases A and B and decide which one is optimal. 

• To measure the uncertainty on the frequency, amplitude, and 
phase. 

The simulated light curves were constructed to have the same 
noise characteristics as for the observations. To measure the proper- 
ties of the noise in the observations, we first subtracted the detected 
frequencies before computing the amplitude spectrum. We used 
two components to describe the white noise and the drift noise: 

(1) The white noise was measured as the mean level in the ob- 
served amplitude spectrum the frequency range 1000-1500 /.iHz. 

(2) The drift noise leads to an increase in the noise level towards 
low frequencies. This "1//" noise was measured by fitting a fidu- 
cial spline function to the observed increase in noise in the ampli- 
tude spectrum. 

Since each site had a different white noise component and a 
different shape of the 1 // noise we constructed the light curve site 
by site before adding the data as we did for the actual observations. 

In practice we created an evenly sampled light curve covering 
the time span of the observed time series, but with an over-sampled 
grid in time by a factor five. Each light curve consist of random 
numbers with different seed number. The mean is zero and the 
numbers have a point-to-point scatter corresponding to the noise 
level found in step (1). Firstly, we take the Fourier transform of 
the light curve to the frequency domain, and multiply the result by 
the fiducial spline fit found in step (2). We then take the Fourier 
transform to return to the time domain, and pick the synthetic data 
points closest in time to the observed data points. Finally, we add 
the frequencies detected in the star. For each star 50 simulations 
were realized. 

We made a Fourier analysis of each simulation as we did for 
the observed data (see Sect. 14. 3| >. In Fig. IB II we show how often 
each inserted frequency was recovered in the simulations versus 
S/N. The top and bottom panels are for EW Cnc and EX Cnc, 
respectively. In each panel open and solid symbols correspond to 
cases A and B (Sect. 14.31 ). For EW Cnc it is seen that in both case A 
and B all frequencies with S/N > 5 are recovered. We prefer case B 
for EW Cnc since it has a significantly lower residual noise level 
(see Table QJ. For EX Cnc we chose case A because it shows a 
much better recovery rate than case B. 

The uncertainties on frequency, amplitude, and phase given in 
Tables [2] and [3] are computed from the rms scatter in the simula- 
tions. We plot the uncertainties versus the amplitude for EW Cnc 
(case B) and EX Cnc (case A) in Fig. lB2l The dashed lines are the- 
oretic al lower estimates of the unc e rtainti es using the expressions 
from [Montgomery & O'Donoghu^ dl999h which are valid for mi- 
correlated noise, ie. for data with only a white noise component. 
For real data one must multiply the un certainty by the square root 
of the estimated correlation length (Montgomer y & O'Don oghue 
1 19991) . We found this by calculating the autocorrelation of the resid- 
ual amplitude spectrum after the mean value has been subtracted. 
The intersection of the autocorrelation function with zero is a mea- 
sure of the correlation length, D. For both EW Cnc and EX Cnc we 
find D = 3.2 ± 0.1 or \f~D = 1.79 ± 0.06. The resulting estimates 
of the uncertainties are marked by the solid lines in Fig. lB2l 

In general we find an acceptable agreement with th e estimated 
uncertainties from [Montgomery & O'Donoghue] dl999l) . For most 
parameters the estimated uncertainties are too low, especially so 
for the frequencies and phases in EX Cnc. This is likely because 
the frequencies lie closer in EX Cnc compared to EW Cnc. The 
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Figure B2. Uncertainties of frequency, amplitude, and phase based on simulations of EW Cnc (leftpanels) and EX Cnc (rightpanels). The dashed lines are 
theoretical lower estimates from Montgomery & O'Donoghue (1999), while the solid lines result from multiplying by the correlation length \/~D. 
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Figure Bl. Recovery rates of frequencies inserted in the simulations of 
EW Cnc (top panel) and EX Cnc. The results for case A and B are shown 
with open and solid symbols, respectively. 



uncertainties on the amplitudes in EW Cnc are close to the estimate 
for uncorrelated noise. 
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